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Abstract 

We consider the three-boson problem with 5-function interactions in one spatial dimension. Three 
different approaches are used to calculate the phase shifts, which we interpret in the context of the 
effective range expansion, for the scattering of one free particle off of a bound pair. We first follow 
a procedure outlined by McGuire in order to obtain an analytic expression for the desired S-matrix 
element. This result is then compared to a variational calculation in the adiabatic hyper spherical 
representation, and to a numerical solution to the momentum space Faddeev equations. We find 
excellent agreement with the exact phase shifts, and comment on some of the important features 
in the scattering and bound-state sectors. In particular, we find that the 1+2 scattering length 
is divergent, marking the presence of a zero-energy resonance which appears as a feature when 
the pair-wise interactions are short-range. Finally, we consider the introduction of a three-body 
interaction, and comment on the cutoff dependence of the coupling. 
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I. INTRODUCTION 



The three-body problem with short-range interactions has been of considerable interest 
for many years in nuclear physics P, S B IE S 0]- More recently, with the realization of 
Bose-Einstein condensates in dilute alkali gases, and the ability to tune the 2-body scattering 
length for such atoms near a Feschbach resonance, it is of increasing interest in atomic 
physics jH 0, EH 11, In light of the relatively recent development of effective field 



theory (EFT), we revisit the old model of three particles in one dimension interacting via 
5-function interactions. While this model has been considered previously by many authors, 
we believe this work provides some unique insights with regard to both atomic and nuclear 
physics. 

Consider the regime where typical nucleon momenta lie well below the pion mass. In 
this limit, it is possible to construct a nonrelativistic EFT in which the pionic degrees of 
freedom are integrated out. This leaves only nucleon fields with contact interactions, and 
higher order derivative corrections. Regarding such an EFT in the two-nucleon sector, a 
great deal of literature has emerged over the past decad e (lil . FLU fl5L 16l | . More recenth 



great deal oi literature nas emerged over tne past decad e 1141 not lipl] . More recently 
there has been a focus on the three-nucleon sector 0, E3, Il8j| . For a recent review see . 
Further, there are now a family of high precision nucleon-nucleon (NN) potential models 
which reproduce NN scattering phase-shifts up to lab energies of 350 MeV (see j2(j and 
references therein). Each of these potentials treat the long-range portion of the interaction in 
the same way via one-pion exchange, but differ in the treatment of the less understood short- 
range physics. Hence, matrix elements of such interactions are said to be model dependent. 
It is possible, however, to decimate the high momentum degrees of freedom by a sequence 
of renormalization group (RG) transformations in order to arrive at a model independent 
low-energy effective interaction j23|. This suggests that low-momentum potential models 
with only nucleons as explicit degrees of freedom may provide a sufficient description of 
few-nucleon systems. Further, EFT may be used to systematize calculations of low-energy 
phenomenon, in principle, allowing calculations of arbitrarily high accuracy. 

For atomic systems, one dimensional Bose gases are of particular interest since phase 
fluctuations are enhanced. It may seem that one dimensional geometries require a radial 
confinement of order the Bohr radius, but this is in fact not the case. All that is required 
is that the energy ga p in the transverse direction be much greater than the gap in the 
longitudinal direction |22j. Also, one dimensional geometries have been observed to display 
higher critical transition temperatures to BEC [23j, and substantially reduced three-body 
recombination rates 24]|. These developments underscore the importance of the three-body 
problem with short-range interactions in one dimension. 

This paper is the first of a pair which investigate EFT and low-momentum effective inter- 
actions in one dimension. For simplicity, we consider only spinless bosons. Scattering theory 
in one dimension plays a central role in all of our calculations. Of particular importance is 
the effective range expansion [i^, |26| , which takes a slightly unfamiliar form. We refer the 
reader to |27[ for the relevant one dimensional derivation. 

We calculate the exact symeterized S-matrix element for the scattering of one boson off 
of a bound pair, and derive an analytic expression yielding the effective range expansion to 
all orders for this 1+2 process. Having found an exact solution, we proceed to calculate the 
adiabatic hyperspherical potential curves in a manner similar to reference |28|. We use the 
eigenchannel R- matrix method [29, 30] in order to determine the scattering phase-shifts, and 
find good agreement with the exact solution. Our results for the phase-shifts, however, differ 



in a critical way from those presented in reference [31| . We trace this disparity to varying 
definitions for the S-matrix element itself. We argue that our definition for the S-matrix 
element is consistent with the threshold behavior of the effective range expansion and with 
the statement of Levinson's theorem in one dimension 3^, 33|. Finally, we derive and solve 



(numerically) the momentum space Faddeev equations for the 1+2 scattering amplitude, and 
find excellent agreement with the exact result. This approach also provides a convenient 
way to analyze the cutoff-dependence of the scattering amplitude and determine the running 
of the three-body coupling constant. 



II. EXACT SOLUTION 



McGuire |24j has shown that when the masses of three identical bosons are the same, 
and the strengths of the pair-wise 5-function interactions are equal, then the elements of 
the scattering matrix can be found by simple geometric optics. Here, we briefly sketch his 
original arguments to calculate the S-matrix, and go one step further to show how this result 
can be examined in the context of the effective range expansion. 

We begin with the interaction expressed in hyperspherical coordinates p and 9 (see ap- 
pendix EJ: 

V(p,9) = {5{\cos9\) + 5{\cos9 - n/3\) + 5{\cos9 + tt/3\)) . (1) 

In the two-dimensional plane covered by p and 9, this interaction is non-zero on three lines 
which intersect at angles of ir/3. Each of the resulting six regions corresponds to a unique 
ordering of the three particles along the real line. The elements of the scattering matrix 
are calculated by tracing a arbitrary ray through the potential diagram, and keeping track 
of the reflection and transmission amplitudes at each intersection. S is then a six by six 
matrix which is indexed by a given ordering of the three particles. The situation is further 
simplified by choosing one particular initial ordering and calculating the six corresponding 
elements indexed by the final ordering. All other elements are readily found by permutations 
of the original ordering. 

We write the familiar transmission and reflection amplitudes as: 

T=-^— (2) 

R=—^- (3) 

with 

2&COS0 . . 

a = = ika 2 cos 0, (4 J 

imco 

where k cos now denotes the momentum component of the initial ray which is normal to 
the surface of the <5-function line. The incoming ray can be traced through the potential 
diagram with the introduction of three angles 0i, 02 and 03 denoting the angle with respect 
to the normal for the first, second and third 5-function line, respectively. If we let Tj = T(aj) 
and Ri = R(ati) be indexed by the wave vector fccos0j, and let the initial ordering of the 
particles be (123) from left to right, then we find the elements of the S-matrix tabulated in 
Table □ 



TABLE I: S-matrix for the 5-function interaction; note there is a common factor in each element, 
and we've defined: 5 = 5- {a\ + l)(a2 + l)(a3 + 1) 



( {outgoing) \S\ {incoming) ) 


Amplitude 


S 


((123)|5|(123)) 


R1R2R3 + T1R2T3 


— 1 — a±a2 


((213)|5|(123)) 


R1R2T3 + T1R2R3 


OL2 


((132)|5|(123)) 


R1T2R3 


0L2 


((231)|5|(123)) 




— Ol\OL2 


((312)|5|(123)) 


R1T2T3 


— OL2OL3 


((321)|5|(123)) 


TiT 2 T 3 


a\Ct20l3 



Boundary conditions for fragmentation states in which two particles are bound by B 2 = 
— are imposed by taking one component of the wave vector to be imaginary, so that a goes 
to —1. For example, if particles 1 and 2 are a bound pair at large p, then we take a\ = — 1, 
so that Ti and R 1 are both divergent. In order to evaluate the scattering amplitude, we 
then are free to set 7\ and R\ to unity while any amplitude not containing either Ti or R\ 
is set to zero. In order to facilitate this, we label the momenta of the individual particles 
with the following kinematics: 

h = i 2 + Te (5) 

h = -\[\q- (7) 

This particular choice satisfies k\ + ki + k 3 = for center-of-mass coordinates, and E = 
^{kf + k\ + fcf) = 2^1 — B 2 , which defines q in terms of the total energy. If we consider 
particle 3 scattering off of a bound state of particles 1 and 2, then there are three available 
options. Either there is total transmission and the ordering goes from [(12)3] to [3(21)] with 
direct amplitude: 

Voga 2 — bt 

or there is rearrangement where the ordering goes from [(12)3] to either [(23)1] or [(13)2], 
each of which occur with the same exchange amplitude: 

A x = T2R3 = ^1 + 4^ (g) 

-3i{qa 2 ) 2 - 4V6(ga 2 ) + 6i 

For the identical boson case, the coherent sum of these three amplitudes yields the desired 
1+2 S-matrix element: 



(10) 



By utilizing the relation: 



we obtain 



,„. n , 2itan5 , . 

exp2z5-l = - — — , 11 

1 — i tan o 

^' = V!(ife (12) 

Expanding this quantity in powers of q 2 yields the effective range expansion to all orders. 
The crucial feature is that the first term equal to the inverse of the scattering length is 
missing, indicating the presence of a zero energy resonance. To be more precise, we would 
expect gtan5 = ^■ + |r 3 (g) 2 + C > ((g) 4 ), but upon inspection of Eq. (JUJ), we see that the three- 
body scattering length is infinite, and the expansion begins with a term 0((qci2) 2 )- It should 
be stressed that this feature persists regardless of the strength of the (^-function interaction. 
There is a state at zero 1+2 collision energy for all attractive zero-range interactions, no 
matter the value of the scattering length. 



III. FADDEEV EQUATION 



The Faddeev approach provides an independent way to analyze the threshold behaviour 
of the scattering amplitude. While many readers are familiar with Faddeev methods, in 
the interest of making the discussion self-contained, we provide a brief derivation of the 
integral equation describing 1+2 scattering in appendix [BJ If we define the amplitude 
K(p, k; E) to satisfy Eq. IB23I with the ie replaced by a principal value prescription, and 



include the normalization of the two-body bound state from Eq. IB101 then we may identify 
kt&nS = —2K(k,k) to obtain an expression which is convenient in the context of the 
effective range expansion. A manifestly three-body interaction parameterized as 2V3/A 2 can 
be included in the kernel in a straightforward manner: 



Z(q,p,E) 



mE — q 2 — p 2 



+ 



2V* 



(mE — q 2 — p 2 ) — p 2 q 2 A 2 



— ~\l —mE + 
3a 2 V 4 



1 + a,2 \ -mE + 



3p 2 



(13) 

In appendix we present an alternative derivation of the kernel above starting from a 
many-body Lagrangian density. 

The numerical solution to principal value integral equations of the form: 



K(p, fc; E) = U(p, k)+V I ^U{p, q)^L_K(q, k; E). 



(14) 



is accomplished by letting q = q n , so that the integral may be written in terms of matrix 
multiplication, and the principal value prescription is enforced by restricting the sum: 



K n ,m U n . m + / dq U m .i o 9 E-ni 



(15) 



Which can be written more succinctly as: 



MK = U 



(16) 
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FIG. 1: The effective range expansion found by numerical solution to the Faddeev equation is 
shown. 



where, 



ml 



"m,l Um,l 



71 ' p 2 n - pf 



K n>i = K(p n ,pi). 



(17) 
(18) 



Inversion of the kernel M is required for each energy (indexed above by n) for calculation 
of the on-shell K-matrix, and hence the phase-shifts. When the interaction has many high 
momentum components, a linear spacing of grid points becomes inefficient, and the inversion 
of the kernel becomes computationally cumbersome. By performing a change of variable 
p = exp(t), it is possible to space the grid points logarithmically, facilitating the solution 
for such interactions. If p = exp (t), then dp = pdt, and the kernel will carry an extra factor 
of the internal momentum q: 



ml 



dt 



Jm,l 



Pi 

ml r, n 

7T Pl - pf 



(19) 



The minimum and maximum momenta may now be chosen to define a domain p e 
[exp (t m i n ),exp (t max )]. The results using the above procedure are plotted in Fig. [T] Clearly, 
as the size of the matrix is increased, the amplitude approaches the exact result of Section ITT1 



IV. ADIABATIC CURVES AND THE EIGENCHANNEL R-MATRIX SOLUTION 



The final apprach involves the adiabatic representation and the eigenchannel R-matrix 
method. We refer the interested reader to appendix El for a review of these tools. For 



5-function interactions, the eigenstates of the adiabatic Hamiltonian inside one of the six 
regions are simply solutions to the free Schrodinger equation. As long as we restrict our 
analysis to one of the six regions, the solution may be found by separation of variables and 
must be of the form l2i 



0o (p, 9) = A (p) cosh (q 6) 
(f) n (p,0) = A n (p) cos (q n 9) 



(20) 
(21) 



We must treat 0o as a special case since it represents the only channel with a two-body 
bound state. The eigenvalue is found by demanding continuity of the wavefunction at the 
5-function surface. If we restrict our analysis to the region 9 e (0, 7r/3), then integration 
of the Schrodinger equation from 7r/6 — e to 7r/6 + e leads to the following transcendental 
equations for q n 



q tanh 



Qn tan 



,g n 7r, 



V2p 
a 2 



(22) 
(23) 



Normalization of <fi n gives: 



Mp) - 



vr+^sinhf^ 

3 . 
7r H sin 



ft 



In 



(24) 
(25) 



The adiabatic potential U n is related to q n by: 



U (p) 

Un(p) 



-<il 



Imp 1 

q 2 
Imp 1 



(26) 
(27) 



We've verified that our potential curves are in agreement with those presented in |2£ 
We provide a plot of the first few in Fig. El As expected, there is only one attractive channel 
which is open below the dimer (two-body bound state) breakup threshold. This model 
supports one true bound state with energy E = —AB 2 . We've found that a calculation with 
one adiabatic channel under-binds this state by about 0.03I?2- Further, the inclusion of more 
coupled adiabatic channels does not serve to improve our result. A calculation without the 
diagonal coupling term gives a second bound state at —1.0021352, very close to threshold. 
When the repulsive diagonal coupling is included, this bound state is no longer supported 
and its eigenenergy is above the two-body binding. The presence of such a state is of course 
consistent with the fact that the exact scattering solution indicates a divergent scattering 
length. A table of these results is shown in Table EJ 

For the case in question, the matrix elements in Eq. IA23I and Eq. IA24I take the form 
(compare with Eq. (|A13|) ): 



A m>n — PoF m (p )F n (p c 



(28) 




FIG. 2: The adiabatic potential curves for the (^-function interaction are shown. 



TABLE II: Eigenenergies for the <5-function model using the adiabatic representation. All energies 
are in units of the two-body binding B2. 



Level 


Exact 


1 channel no Qoo 


1 channel with Qoo 


5 channels full calculation 


Faddeev 


1 


-4 


-3.96902 


-3.96106 


-3.96106 


-3.9998 


2 


-1 


-1.00213 


-0.99587 


-0.99587 


-1.0000 



Po 

,„.„ - , pdp S, 
'0 

"Po 

pdp 

'0 



OF dF 
dp dp 



+ F m F n {k 2 - 2mU n {p)) 



,91,,, 

Pm,n\P)Fm I 2— -\- J F n -\- QmnF m F n 

op p 



(29) 



We choose a set of 7 th order b-splines as our basis set in the expansion for the the functions 
F n (p). Basis splines have prov en to be a versatile and efficient basis set for a wide variety 



systems |35|, |3y, |37 



See |38( for mathematical details, fast algorithms and fortran code. 
The results presented in Fig. |3] and Table |Il] are for a set of 40 splines with a quadratic 
distribution of knot points over the region p e [0, 20 02]. 

In the scattering sector, there are a number of subtleties involved with the matching of 
the wavefunction in the asymptotic region. If the diagonal coupling term is ignored, then 
the Schrodinger equation in the ground state channel for p — > 00 takes the form: 



2 



1 d 2 



dp 2 p dp 2 



F (p) = 



(30) 



This is recognized as the zeroth order Bessel's equation, and hence / must be of the form: 



F (p) = AJ (qp) + BY (qp) (31) 
where q 2 = k 2 — k%. If the diagonal coupling is included, then the solution must satisfy 

d 2 1 d 2 



dp 2 p dp 2 



+ k 2 - k 2 - Qoo(p) 



F (p)=0. (32) 



Qoo behaves asymptotically as l/4p 2 , meaning that the solutions must now be fractional 
order Bessel functions, 

F (p) = AJ 1/2 {qp) + BY 1/2 {qp). (33) 
Reference j3l[ defines the scattering matrix in terms of the outgoing wavefunction as 

Fo(p) - VW)(H[%(qp) + e™'H[f 2 (qp)), (34) 
which, ignoring overall normalization, and taking qp 3> 1, can be written 

■F (p) -> e~ iqp - e 2i5 'e iqp (35) 



Reference |3jj asserts that the minus sign appearing above adds an extra it/ 2 to the phase- 
shift, so that the total phase shift starts at 37r/2. This assertion would indeed be consistent 
with Levinson's theorem in three-dimensions, where one would obtain a ir from the known 
bound state and a 7r/2 from the zero-energy resonance. However, in one spatial dimension, 
we note that there is no additional tt/2 for the zero-energy resonance. The statement of 
Levinson's theorem in one dimension for the even parity solution takes the form 



lim5 e = (n e — l/2)7r noncritical case (36) 

lim S e = n e Ti critical case. (37) 

The critical case applies when there is a zero-energy resonance. The statement for the 
odd-parity solution is identical to that for three-dimensional S-waves: 

lim5 = n 7r noncritical case (38) 

lim5 = (n + l/2)n critical case (39) 

We are concerned only with the even parity case. We define our scattering matrix in the 
following fashion. We require S to be the coefficient multiplying the outgoing e iqp in the 
limit qp ^> 1. In terms of the out-going wavefunction, this means 

F (p) -> (e- iqp + e 2iS e iqp ) . (40) 



We justify our choice by considering the limit £i <C £2; which is appropriate when particles 
1 and 2 are bound and particle 3 is far away, in which case, p~£2(l + !%"---)- m this way, 
the product F o (p)(j) o (p,0) represents a 2-particle bound state in one relative coordinate and 
an oscillatory wave in the second relative coordinate, symeterized over all permutations of 
particles. Note also that the extra factor of \j-Jqp appearing in our asymptotic solution 



is different from the convention of reference |3lj : this is because we choose to work with 
the full wavefunction instead of the reduced wavefunction, hence our integration measure 
remains J pdpdO. 

In terms of standing wave solutions, the definition of the phase-shift above is consistent 
with the following expression involving tan (5): 

Fo(p) — > ( cos (qp) — tan (8) singp) (41) 



The form in Eq. (f4*T]) leads to the desired product state corresponding to 1+2 scattering in 
one dimension. By identifying 

b = _d)nF (p) 
dp 

one may easily solve for tan 5. An alternative argument may be formulated by simply noting 
that the normalization condition J |^0o(p, 6) = 1 requires that 0o scale like y/p. This means 
that the full wavefunction is proportional to the quantity ^/pF (p). If we consider the 
wavefunction at a particular angle, and let ^(p,6 = 6o) — > Acos (qp + 5), which is the 
proper asymptotic form for an even solution, then 

d\n^(p,e = e Q ) ,1 . 

^ = - h + Y P = ~ q tan (gp + } ' (43) 

This expression is entirely equivalent to the matching condition Eq. 021 We note that our 
phase shift 5 is related to 5' by 

tan5=^ 1 — (44) 

tan 5' v 1 

Clearly, this will alter the behavior at threshold. We note that with our definition, the 

presence of a zero energy resonance is consistent with the threshold behavior of the effective 

range expansion, namely that the 1+2 scattering length is divergent. 



V. DISCUSSION 

The most striking feature of these results is the presence of the zero-energy resonance 
marked by the divergent scattering length. As the cutoff is lowered and the range of the in- 
teraction becomes finite, the scattering amplitude in the limit k — > becomes nonzero. This 
behavior is illustrated in the off-shell amplitude calculated numerically via Eq. IB23I (except 
with a principal value prescription) shown in Fig.QJ which is the ID analog of Fig. 5 appear- 
ing in reference Q. There is clearly a fixed point in the limit A — ► oo. The presence of the 
zero-energy resonance is further substantiated by the variational calculations of Section II VI 
When the repulsive second- derivative coupling Qo,o(p) is omitted, the interaction supports a 
second bound state with energy E = — 1.00213I?2- When Qo,o(p) is included, the bound state 
is no longer supported and a solution to Eq. IA12I gives E > —B 2 , consistent with the upper 
bound theorem. These results are also consistent with the threshold behaviour of qa 2 tan S 
shown in Fig. El The Eigenchannel R-Matrix method gives lim ga2 _,.o QO>2 tan S ~ 0.065 when 
the diagonal coupling is excluded, and lim ga2 ^ Q a 2 tan 5 ~ —0.008 when it is included. 
Finally, the exact effective range expansion calculated in Section |H] indicates a divergent 
scattering length in perfect agreement with the numerical calculations. It is again of consid- 
erable note that the zero energy resonance is present regardless of the value of the 5-function 
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FIG. 3: The effective range expansion for 1+2 scattering with <5=function interactions calculated 
using the Eigenchannel R-Matrix method is shown. A calculation with five coupled channels is 
well converged such that the inclusion of more channels does alter the results. 



coupling, or equivalently the two-body scattering length. It appears at exactly zero relative 
energy in the 1+2 system as long as the two-body interactions are of zero range, and moves 
away from threshold as the interactions become of finite range. 

It is now instructive to consider the on-shell results when various cutoffs are enforced in 
Eq. IB23I As the cutoff is lowered, we require that 1/03 = (03 being the 1+2 scattering 
length). This quantity is proportional to the scattering amplitude satisfied by the principal 
value version of Eq. IB231 With the introduction of a three-body interaction, the cutoff 
dependence of the amplitude can be absorbed into the three-body coupling V3, yielding a 
largely cutoff invariant amplitude as seen if Fig. 0] and Fig. 

This paper has treated three bosons only at zeroth order in EFT. In a second paper, 
we shall extend the analysis to include the two-body effective range and shape parameter. 
Predictions from the resulting EFT will be compared to calculations using a realistic phe- 
nomenological NN interaction. Three-nucleon observables will serve as a testing ground for 
the effective theory. 



Acknowledgments 

NPM would like to thank Chris Greene for many helpful discussions regarding the eigen- 
channel R-matrix method and the adiabatic hyperspherical representation, Daniele Bor- 
tolotti for discussions on b-splines, and the INT in Seattle for their hospitality. NPM and 
JRS would like to thank Jae Park for discussions on ID Bose gases. 



T 




FIG. 4: The off-shell amplitude for our one dimensional model is shown. The black curve which 
should be considered the "exact" result is for Aa2 = 10 4 . As the cutoff is lowered to Aa2 = 5, a 
three-body term with V3 = —0.0516 brings K(0, 0) into agreement with the exact result. As the 
cutoff is lowered further to Aa2 = 2, a three-body term with V3 = —0.17595 is required. Finally, 
V3 = —1.907 produces the exact K(0, 0) for Aa2 = 1. 



APPENDIX A: REVIEW OF THE ADIABATIC HYPERSPHERICAL REPRE- 
SENTATION AND THE EIGENCHANNEL R-MATRIX METHOD 



The essential strategy for solving the Schrodinger equation in coordinate space is to 
transform the partial differential equation (PDE) into a set of coupled ordinary differential 
equations (ODEs). While there are a variety of representations that realize this goal, the 
one best suited to the present problem is the adiabatic hyperspherical representation [3^ 
(see also (HI). 

Let us first introduce the appropriate relative and hyperspherical coordinates. Jacobi 
coordinates in one dimension for equal masses are denned via 




/ _L JL J_ \ 

1 73 73 73 » 

J- =1 

72 72 

\ 76 76 76 / 




(Al) 



where X marks the position of the total center-of-mass, £1 is the relative coordinate for the 
first two particles , and £2 is the relative coordinate between third particle to the center-of- 
mass of the first two. Hyperspherical coordinates in one dimension are simply circular polar 




(qa 2 ) 2 

FIG. 5: The on-shell amplitude for our one dimensional model is shown for various values of the 
cutoff. As the cutoff is lowered, a three-body contact interaction of natural size is introduced to 
fix the 1+2 scattering length. 




coordinates: 



£ 1 = pcos# (A2) 
£ 2 = psin0 (A3) 

p 2 = & + e 2 (A4) 

Here, p is a measure of the general size of the system. At small values of p, all three particles 
are in close proximity, while at large values of p, the situation depends on the angle 9. There 
are some values of 9 that correspond to two of the three particles being near each other, and 
other values of 9 where all three particles are far apart. 

With this transformation, the Schrodinger equation becomes (compare with Eq. 

li{%) + 7^ +k2 - u(p ' e) h p - e) =° (A5) 

where k 2 = 2mE and 

U(p,0) = 2m(v(V2p\cos(9)\) + 

V(V2p\ cos(9 + 7r/3)|) + V(V2p\ cos(9 - tt/3)|)V (A6) 

This potential has a very high degree of symmetry: 

U(p, 9) = U(p, 0±\) exchange (A7) 

U(p, 9) = U(p, -9) parity (A8) 

The combination of exchange and parity result in a six-fold symmetry allowing the angular 
part of the wavefunction to be represented as a sum over terms proportional to cos 6n9. 

The reduction of the PDE Eq. IA5I into a set of coupled ODEs is accomplished by expand- 
ing the wave function into a sum over different adiabatic channels: 

*(p,9) = J>nW) = J2Mp,0)FM- (A9) 

n n 

where n are defined as eigenstates of the adiabatic Hamiltonian: 

H ad (p, 9)4> n (P, 0) = k 2 n (p)<f> n (P, (A10) 

with 

H ad (p,0) = ^r-^ + U(p,9). (All) 

It is important to note that Eq. lAlOl depends only parametrically on p. We impose boundary 
conditions such that n is even at 9 = and 9 = tt/6. Inserting the expansion in Eq. (|A9|) 
into a variational expression of the form 

_ J v dvn-v 2 + u)* (AU) 

k ~ JydVW (A12) 



and demanding that the solution be stationary with respect to variations in the functions 
F^, one arrives at the following matrix equation: 



(T + U)F = k 2 F. 



(A13) 



We've defined 



d 2 2 d \ f d 2 s 

dp 2 pdpj V dp p/ 



and 



U m , n (p) = 5 m , n {k 2 m (p)). 
Also, we've introduced the non-adiabatic channel couplings defined as 

d , 



Qm,n{p) = U>r, 



dp 



d 2 



dp 2 



(A14) 
(A15) 

(A16) 
(A17) 



The solution of the the adiabatic equatation Eq. IA10I Eq. (jAllj) and the evaluation of the 
first and second derivative couplings accounts for the vast majority of the computational 
effort in solving the three-body problem using this approach. The first- derivative couplings 
Pm,n can be evaluated using a Feynman-Hellman like argument. For a p-parameterized 
system defined by H(p)<f) = E(p)<p, the Feynman-Hellman theorem states 



dE 
dp 



OH 



dp 



(A18) 



This relation can be used to find that 



d_ 

dp 



(<t>m If 0n> 



dp 



k n(p) ~ k li(p) 



(A19) 



for m ^ n, and vanishes for m=n. The second-derivative couplings, Q m ,ni may be readily 
calculated by noting: 



Qm,n \P \m,n 



dP 
dp 



(A20) 



Use of the above relation, however requires calculating the first-derivative couplings between 
many channels so that the square of P converges. For calculations involving a single adi- 
abatic channel, it is more convenient to estimate Qo,o by using a 3-point or 5-point rule. 
This involves solving the adiabatic Hamiltonian at three consecutive nearby values of the 
parameter p, and calculating the second derivative numerically before evaluating the inner 
product. 

In order to calculate wavefunctions and scattering amplitudes in the scattering sector, 
we use the Eigenchannel R-matrix approach [29, 30], which is a variational calculation for 
minus the log-derivative of the wave function on the surface S of some reaction volume V. 
More precisely, this method finds variational solutions that have a constant log- derivative 
on the surface such that ^ + = 0. 



Starting with Eq. IA12| we define b = — 9lr ^^ , where n represents the unit normal vector 
to the reaction surface S (h = p in our case); application of Green's theorem to the kinetic 
energy term in Eq. (jA12|) allows us to write an expression for b at a fixed k 2 . 



f dV [(-V*t • V*) + ^(k 2 - 17)* 
b= l — 



J dS M 



(A21) 



Note that we were able to factor b out of the surface integral in the denominator only because 
the desired solution has a constant log-derivative on the surface 5*. Eq. (jA21|) is an identity 
obeyed by exact eigenstates of the Schrodinger equation Eq. (jA5|) that have a constant b on 
S. By taking the first-order variation of this expression with respect to small deviations in 
this expression can be shown to be a variational expression for b. 

In the adiabatic representation, we expanded the wavefunction according to Eq. IA9I Now 
we expand F n = c nt0l B a , and Eq. (jA21|) is cast into the form of a generalized eigenvalue 
equation: 

bAc = Tc 6^ - ' (A22) 

op 

where 



A m ,„ = / dS -n ip m *Pn (A23) 
Js 

T m ,n= [ dV[-V4> m ■ V^n + *Pm(k 2 - U)*P n ] (A24) 

Jv 

APPENDIX B: DERIVATION OF FADDEEV EQUATIONS 

In this section, we derive the properly symeterized integral equation satisfied by the 
scattering amplitude for one free boson off of a bound pair. Our derivation relies heavily on 
the original work of Faddeev Lovelace and Amado jUj], however we largely hold 
to the notation conventions of Watson and Nuttall 44]| . We shall begin by calculating the 



two-boy scattering amplitude for a separable interaction of the form: 

V = c \g)(g\, (Bl) 

with momentum space matrix elements: 

W\V\q) = Cog{q')g{q). (B2) 

The T-matrix element (q'\T(E)\q) satisfies the Lippman-Schwinger equation: 

rp/ i ^ ( a / x , / / \ f dq" mg{q")T(q",q,E) 

T(q ,q,E) = c Q g{q )g(q) + c g(q ) J — m ^ _ qll2 (B3) 

The solution is found by defining the energy dependent function 3 



which satisfies the algebraic equation: 



h(q,E) 



dq" mg 2 (q") 
2tt mE - q" 2 



[cog(q) + c Q h(q,E)}. 



(B5) 



Solving for h(q, E) and substituting the result into the Lippman-Schwinger equation quickly 
yields the solution: 



T(q',q,E) = g(q')g(q) 



1 

Co 



dq" mg 2 {q") 
2tt mE - q" 2 



(B6) 



Taking the limit g(q) — > 1 is equivalent to solving the corresponding Schrodinger equation 
with contact interactions. This yields a T-matrix element independent of q and q': 



mT{E + ie) = mco 



1 + 



mc 



1\l — mE — ie_ 



(B7) 



In the center of mass frame, the energy is written E = k 2 /m and the scattering amplitude 
/ for the even parity wave is related to the on-shell T-matrix by: 



/ = £ (*» ~ 1) 



itanS mT(k,k,E) 



2ik 



1 — it&nS 



2ik 



1 _ mc 
2ik 



(B8) 



The cross-section is a normalized probability in one dimension and is given by a = \f\ 2 . 

If the coupling is negative, the interaction supports a bound state. The T-matrix element 
will exhibit a pole at the binding energy E = —B 2 . From inspection of Eq. ()B7|1 . it is clear 
that this requires B 2 = mc^/A = 1/ma 2 . The state vector for the bound state is: 

\iP B ) = VNG (-B 2 )\g). (B9) 

where G (—B 2 ) = (—B 2 — H )~ 1 is the free particle propagator evaluated at the bound state 
energy. The normalization constant is easily evaluated by contour integration to be 



N 



-i 



( 9 \Gl(-B 2 )\g) 



2-3 



ma 



(BIO) 



The indices in the three-body sector follow the convenient "odd-man-out" notation. 
When considering matrix elements of two-body operators in the three-body state-space, 
the two-body operator V a with a = 1 shall denote the interaction between particles 2 and 3. 
We are only concerned with internal degrees of freedom, and will therefore work in the total 
center of momentum frame pi + p 2 + p% = 0. Matrix elements will be taken with respect 
to state vectors of the form \p a ,q a ), where p a represents the momentum of the spectator 
particle a, and q a represents the relative momentum of the remaining two particles. Let 
|V£> = + |4 2) ) + |4 3) ) describe an eigenstate of the full Hamiltonian H which corre- 
sponds to an initial state \\ a ) with the two particles not equal to a forming a bound state. 
The solution is found by solving the Faddeev equations for the components 



|Vf ) = G (E)\a, Pa )8 aP + G (E)T^}. 



(BH) 



We write the two-body T-matrix in the three-body state space as: 

,2 



T ? {E) = j if l/? ' P " >7 * { E ~ ^ (B12) 

where Tp(E) is the dimer propagator in the f3 channel: 



rp{E) 



'_! /" dq rng}{q) 



Cq J 2ir q 2 — mE 



(B13) 



For gp(q) — > 1, Tp(E) is equal to the two-body T-matrix found in the previous section, left 
multiplying Eq. (|B11J) by (<5,p<s| and summing over [3^5 leads to: 

X S)a (ps,Pa) = Z 5) a(ps,p a ) -^ Z sAV5,Vi3)rp \ E - -£^J Xp >a (pp,p a ). (B14) 

The amplitudes X^tj^, p Q ) and Z^ a (ps,p a ) are defined as: 

X 5ia (p 5 ,p a ) = ^(5,p 5 |Vi 7) ). (B15) 

ZsAPs,Pa) = (1 - <^ Q )(<5,p 5 |G (-£)|a,Pa)- (B16) 

The Born amplitude Z$ tCl (ps,p a ) describes the interaction mediated by the exchange of a 
single particle, and requires calculating q$ in terms of p a with 5 ^ a. The kinematics for a 
given case must be determined by cyclic permutation of the particles j42|. For example, 



^2,l(P2,Pl) " 3p? 1 I 2 • (B17) 

For identical bosons the quantity of interest is the symeterized amplitude given by the sum of 
the direct and exchange pieces X(p, k; E) = X^ D \p, k; E) + 2X^ N \p, k; E). This amplitude 
satisfies the following integral equation: 

X(p, k; E) = 2Z(p, k; E) + 2 J ^-Z(p, q; E)r (e - X(q, k; E) (B18) 

where the Born term for g(q) — » 1 is given by: 
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Z(p, q; E) = — - - , (B19) 

mE — q z — p z — q ■ p 

and the total energy is E = ^-—B 2 . Next we perform an angle average over the dot product. 
In one dimension the angle average of an arbitrary function f(p-q) is / = ^(f(jpq) + f(—pq)), 
and so the Born amplitude becomes: 

rp _ 2 _ 2 

z (p> ?; E ) = -r-p — 2 2 / 2 2 - ( B2 °) 

[mE — q z — p z ) z — p z q z 



The desired amplitude now satisfies the integral equation: 
X{p, k,E + ie) 



2m 



Z(p,h,E) 



dq 



+ 4 / ^-Z(p,q,E) 



-i -i 



-a 2 + 



-mE + - ie 



X(q,k,E + ie) 
2m 



• (B21) 



It is desirable to remove the pole in the dimer propagator and bring this equation into the 
form of the Lippman-Schwinger equation; to this end we define the amplitude: 

-i -i 



X(p,k,E) 
k 2 — p 2 



-a 2 + 



1 



X{p,k,E) 
2m 



(B22) 



A bit of algebra shows that the new amplitude satisfies the equation: 

X(p, k,E + ie) = Z{p, k, E + ie) — — [°° dq Z(q, p, E + ie) X ^ k,E + ^ 
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q 2 — k 2 — ie 



with 

Z(q,p,E) 



mE — q 2 — p 2 



(mE — q 2 — p 2 ) — p 2 q 2 



4 



3p 2 



— rW —mE H 

3a 2 2 V 4 



1 + a 2 



-mE H 

4 



(B23) 



(B24) 



It is computationally more convenient to deal with an amplitude which is real below the 
breakup threshold by writing the above integral equation in terms of a principal value 
prescription using the well known formula: 

= Vt in6(w) (B25) 

uj ± %e 



APPENDIX C: ALTERNATIVE DERIVATION OF EQ. IB181 

Bedaque et al 0] have considered the 3-D three-body problem with short range inter- 
actions in a ground-breaking paper. It is instructive to consider their approach in a 1-D 
context, and that is the purpose of this appendix; a complete analytic sum of the series 
arising in perturbation theory has been found by Thacker 4^ . 



Consider the Feynman rules resulting from the Lagrangian density: 

C = <p(x) (id + 4>{x) - | {^{x^ix)) 2 - ^ (0 t (^)0(^)) 3 • (CI) 

Let us first sum the perturbative series of bubble diagrams for the two-body problem with 
do = 0. Let (po,p) = (^(ki — k 2 ) 2 , ki + k 2 ) with \k\\ = k denote the two- vector in the center 
of momentum frame. The following loop integral is readily evaluated by contour integration: 

L = I ^ * * (C2) 

_ -im _ m 

^=T " 2l' {L6) 



2 x /-mpo + \ - ie 



FIG. 7: The dressed dimer propagator determined by the geometric series of loop insertions. 




FIG. 8: The diagrammatic representation of the integral equation for the 1+2 scattering amplitude 
is shown. 



and the geometric series is easily summed to reproduce the result of Section [Bj 

iA = -iT = -ic (1 + (-ic L) + Hc L) 2 ...) = = (C4) 

Kaplan suggested that the Lagrangian Eq. IC1I may be conveniently rewritten in terms 
of a dummy field D: 

C = 0t (id + ^] + AD ] D - 4=P f # + <fi$D) + h(D ] D<^<p + <jxpDD*) (C5) 
V 2m/ V2 

Gaussian path integration over the auxiliary field D shows that the couplings appearing in 
Eq. ((OSJ) are related to those in Eq. (JUT} by g 2 / A = c and -3hg 2 /A 2 = d . The bare 
dimer propagator is if A, while the sum of diagrams shown in Fig. [7| yields 

,A(^^(l + (^) + (^) 2 + ...) (C6) 

(C7) 



A+™f(- mpo + £l-ie)-V2 

For the 1+2 integral equation, we choose the same kinematics as Bedaque et al Q. 
Let the incoming particle and dimer have two- momenta (k 2 /2m, —k) and (k 2 /4m — B 2 ,k), 
respectively. The outgoing particle and dimer are off-shell with two-momenta (k 2 /2m—e, —jp) 
and (k 2 /Am — B 2 + s,p), respectively. Our integral equation is identical to Eq. (5) in [7[, 
except that the integration measure is J for 1+1 dimensions: 

it(k, p, e) = - 2g 2 iS(-k 2 /4m - B 2 + e,p + k) + ih 

p-^S(k 2 /2m ~ £ ~io, -q) 

[(-2g 2 )iS(-k 2 /4m - B 2 + 2e + q ,p + q) + ih] 
iA(k 2 /4m - B 2 + e + q , q)it(k, q,e + q ) (C8) 



The J ^ integral is readily evaluated by contour integration since the two nucleon propa- 
gators have poles in opposite half-planes. The result is: 



it(k,p) —%2mg 



+ 



h 



— i2mg / — 
2ix 



k 2 + p 2 — mE + p ■ k — ie 2mg 
dq 



+ 



h 



q 2 + p 2 — mE + p- q — ie 2mg 2 
t{k,q) 



A + ^fff -mE-it}- 1 / 2 



(C9) 



With the chosen kinematics, the total energy is E = 3k 2 /4m-B 2 . We've set e = and 



as in reference defined t(k : p) = t(k,p, 2 P )■ Now, averaging over the p ■ k brings us to: 

k 2 + p 2 — mE 



t{k } p) =2mg z 



(k 2 + p 2 — mE) 2 — p 2 k 2 
q 2 + p 2 — mE 



+ h 



2ir [ (q 2 + p 2 — mE) 2 — p 2 q 2 

tj&q) 

a 2 - [3/4{q 2 - k 2 ) + mB 2 - ie}- 1 / 2 



+ 



2mg 2 



(CIO) 



The on-shell amplitude must include the wave-function normalization for the two-body 
bound state, which in field theory is conventionally written: 



with 



dpo 



mg 



-m 



2 21 I 

m g \a 2 \ 



(Cll) 



(C12) 



po=-b 2 2 2 (mB 2 )-y 2 4 

It is desirable to bring this equation into the form of the standing-wave Lippman- 
Schwinger equation. To this end, we define the function a(k,p): 



a(k,p) 



t(k,p)/2mg 2 



p- 



k 2 a 2 - (3p 2 /4 - mE)- 1 / 2 



(C13) 



Since the integral is even in q (indeed, it is only a function of q 2 ), the limits may be taken 
from zero to A provided that we multiply by an overall factor of 2. This of course introduces a 
sharp cutoff A. The integral equation is now written in terms of a principal value prescription 
as: 



4 

a(k,p) = M(k,p; E) V 
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where the kernel M(q,p; E) is defined: 

mE — q 2 — p 2 



A dQM(q,p;E)^ 



(C14) 



M(q,p; E) 



(mE — q 2 — p 2 ) 2 — q 2 p 2 



2mg 2 



4 

3an 



y/-mE + 3p 2 /A (l + a 2 ^/-mE + 3p 2 /4^ 



(C15) 



It is now clear that a(k,p) satisfies the same integral equation as K(k,p;E). 
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